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Simple  models  play  a key  role  in  the 
microstructural  analysis  of  mechanical 
failure  in  composites  and  other  materi- 
als with  complex  and  often  disordered 
microstructures.  Although  equal  load- 
sharing-models are  amenable  to  rigor- 
ous statistical  analysis,  problems  with 
local  load  enhancements  near  failed  re- 
gions of  the  material  have  so  far  re- 
sisted exact  analysis.  Here  we  show  for 
the  first  time,  that  some  of  the  simpler 
of  these  local-load-sharing  models  can 
be  solved  exactly  using  a sub-stochastic 
matrix  method.  For  diluted  fiber  bun- 
dles with  local  load  sharing,  it  is  possi- 
ble to  find  a compact  expression  for 
the  eharacteristie  equation  of  the  sub- 
stochastic matrix,  and  from  it  obtain  an 
asymptotic  expansion  for  the  largest  ei- 
genvalue of  the  matrix.  This  in  turn 
gives  the  asymptotic  behavior  of  the 


size  effect  and  statistics  of  the  fiber- 
bundle  models.  We  summarize  these 
results,  and  show  that  the  important 
features  of  the  exact  result  can  be  ob- 
tained from  a single  scaling  analysis  we 
had  developed  previously.  We  also 
compare  the  statistics  of  fracture  with 
the  appropriate  limiting  extreme-value 
survival  distribution  (a  Gumbel  distribu- 
tion), and,  as  previously  indieated  by 
Harlow  and  Phoenix,  note  that  the 
Gumbel  distribution  performs  quite 
poorly  in  this  problem.  We  comment 
on  the  physical  origin  of  this  dis- 
crepancy. 
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1.  Introduction 

It  has  been  known  since  the  pioneering  work  of 
Chaplin  [1]  and  well  known  since  the  classic  work 
of  Griffith  [2]  that  randomly  occurring  flaws  or 
weak  links  effectively  determine  the  observed  ten- 
sile strength  of  materials.  Early  on  it  was  realized 
that  the  dependence  of  failure  upon  the  weakest 
part  of  a material  structure  gives  rise  to  non-Gaus- 
sian  statistics  for  fracture  stress  and  strain.  These 
developments  lead  to  the  classical  period  of  the  de- 
velopment of  the  statistics  of  extremes  by  mathe- 
maticians such  as  Dodd  [3],  Frechet  [4],  Fisher  and 
Tippett  [5],  von  Mises  [6],  Gnedenko  [7],  and 
Gumbel  [8]. 


Following  the  work  of  Duxbury  et  al.  [9-11], 
there  have  been  many  attempts  to  use  random  net- 
work models  to  determine  the  statistics  and  size 
dependence  of  material  breakdown  [12-16].  These 
calculations  have  in  many  cases  elucidated  the  gen- 
eral behavior  and  size  dependence  of  breakdown, 
but  few  exact  results  have  been  produced. 

Perhaps  the  simplest  model  that  shows  the  statis- 
tics of  brittle  failure  has  been  the  pure-flaw,  chain- 
of-bundles  model  of  Harlow  and  Phoenix  [17] 
which  has  been  studied  by  Harlow  [18]  and  more 
recently  by  Harlow  and  Phoenix  [19].  In  this  model 
there  is  a series  or  chain  of  m structurally  and 
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statistically  independent  bundles  of  n elements 
each  as  shown  in  Fig.  1,  where  the  vertically  ap- 
plied uniaxial  stress  is  shared  by  the  surviving  verti- 
cal fibers  (bonds).  Each  element  or  fiber  is 
independently  present  with  probability  p and  ab- 
sent with  probability /=  1 —/?.  The  survival  proba- 
bility of  the  chain  of  bundles  is  then  the  survival 
probability  of  a single  bundle  raised  to  the  power 
m.  The  main  difficulty  in  this  analysis  is  calculating 
the  survival  probability  of  a single  bundle.  The  ex- 
tension of  this  theory  to  the  survival  of  two-dimen- 
sional networks  is  straightforward  and  amounts  to 
the  approximation  that  cracks  or  flaws  only  exist 
and  break  along  the  direction  transverse  to  the  di- 
rection of  the  applied  stress.  Following  Harlow  and 
Phoenix,  we  assume  the  local-load-sharing  model 
for  a flaw  (crack)  of  length  n to  be 

cr,ip  = o-(,(H-n/2),  (1) 

which  is  to  say  that  the  entire  force  applied  to  the 
cluster  is  concentrated  at  the  tip  (on  the  fibers  ad- 
jacent to  each  end  of  the  flaw  or  vacancy  cluster). 
Failure  of  any  surviving  bond  (all  of  which  have  the 
same  strength)  leads  to  a rip  which  causes  failure 
of  the  entire  bundle.  Solution  of  this  model  re- 
quires finding  the  bond  (weakest  link)  which  expe- 
riences the  largest  stress  enhancement  and  that 
stress  which  would  break  this  most  stressed  bond. 


• • • • • • • 


Fig.  1.  A one  dimensional  array  of  intact  bonds  (fibers)  and 
flaws  (vacancies).  The  tensile  stress  a,  is  applied  vertically. 

Duxbury,  Leath,  and  Beale  [11]  showed  how  a 
one-dimensional  model  such  as  this  could  be  used 
as  a simple  model  for  fracture  or  breakdown  of  a 
two-dimensional  network.  If  one  considers  that 
cracks  or  flaws  only  exist  and  break  horizontally, 
then  the  two-dimensional  model  becomes  that  il- 
lustrated in  Fig.  2 (i.e.,  no  horizontal  bonds  break). 
Then  if  we  impose  spiral  boundary  conditions, 
where  the  last  site  in  a row  is  connected  to  the  first 
site  of  the  previous  row  on  the  other  side  of  the 
sample,  then  the  NxN  network  problem  is  re- 
duced to  a one-dimensional  chain  of  fibers  (or 
bonds)  in  parallel  like  that  in  Fig.  1. 


Fig.  2.  A two  dimensional  lattice  with  only  horizontal  cracks. 
Spiral  boundary  conditions  identify  each  site  on  the  right  edge 
with  the  site  on  the  left  edge  of  the  previous  row. 

2.  Single-Cluster  Calculation:  The  Sub- 
Stochastic  Transition-Matrix  Method 

In  this  calculation  one  assumes  that  the  weakest 
link  is  the  pure  flaw  or  cluster  of  vacancies  of  the 
largest  size  that  exists  in  the  sample.  The  survival 
probability  is  then  closely  related  to  the  probability 
that  in  a sample  of  length  L(—N^)  that  there  are 
no  clusters  of  vacant  bonds  (flaws)  of  size  greater 
than  some  prescribed  value  n . Using  a generating- 
function  technique  Duxbury,  Leath,  and  Beale  [11] 
calculated  exactly  the  asymptotic  form  of  the  prob- 
ability to  be 

CL{n)^[\-pr'f  (2) 

in  the  limit  of  large  L . It  is  now  possible  to  rederive 
this  result  while  introducing  the  sub-stochastic 
transition-matrix  method.  Following  Harlow  [18] 
we  define  all  possible  endings  of  a fiber  bundle  of 
length  L + \ and  the  way  in  which  those  endings 
may  be  generated  from  a bundle  of  length  L and 
the  probabilities  of  those  endings.  Since  there  are 
no  allowed  clusters  larger  than  n,  the  allowed  bun- 
dle endings  or  distinctive  endings  at  a particular 
site  are  an  occupied  site  (1)  followed  by 
vacancies  (0)  so  that  these  distinct  endings 
are  spanned  by  the  basis  vector  <j)L^  = 

(</)(i),</>(io),</)(i(x))....,</>(i()..()))l  where  the  last  element 
contains  n zeros. 
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Then  the  probability  increment  for  going  from  a 
cluster  of  size  r to  a cluster  of  size  r'  on  the  next 
site  is  included  in  the  matrix  product 


p p p ...  p 

</’(!) 

<A(1) 

/ 0 0 ...  0 

</>(10) 

<b(W) 

0/  ••.  ■•.  i 

= 

: : 0 0 

0 0 ...  / 0 

<A(I0...0) 

/ 

which  is  the  same  as  the  matrbc  M operating;  times 
on  the  probability  vector  0o  for  the  starting  site,  or 

(3b) 

The  probability  Ciin)  that  there  are  only  clusters 
up  to  size  n in  the  entire  bundle  (or  network)  of 
size  L is  thus 

CL{n)  = X{4>e)L,  (4) 

e 

where  the  sum  is  over  all  the  elements  of  <j>L.  The 
simplest  and  most  natural  boundary  conditions  are 
periodic  ones  where  Cz.(n)  becomes  the  trace 

CLin)  = tr(M^)^XXt,  (5) 

k 


since  the  1st  and  Lth  sites  must  be  the  same,  where 
\k  are  the  eigenvalues  of  M.  We  find  the  eigenval- 
ues of  M via  its  characteristic  equation 


(u-1) 

a 

...  a 

a 

b 

-1 

0 ... 

0 

D„=detlM/A-/|  = 

0 

b 

-1 

0 

= 0, 

0 

0 b 

-1 

(6) 

where  a =pl\  and  b =//A 

= (1 

-p)/\. 

A cofactor 

expansion  of  the  determinant  D„  about  its  last  row, 
gives  immediately  the  recursion  relation 

= -D„-i  + (-l)''ufo",  (7a) 

where  D„-i  is  the  n xn  determinant  for  clusters  up 
to  size  {n  - 1).  With  Do  = (a  - 1),  the  solution,  upon 
iteration  of  Eq.  (7a),  is  the  characteristic  equation 

{-ly  Dr,={a  -l)  + ab  +ab^  + - + ab"  = 0.  (7b) 

Summing  this  geometric  series  we  obtain 

+ + (8) 

This  equation  is  the  characteristic  equation  of  M 


times  (A  — /),  so  there  is  an  additional  spurious 
root  at/ (since  we  are  interested  in  the  largest  root, 
this  does  not  affect  the  analysis).  Since  M is  primi- 
tive and  non-negative,  its  largest  eigenvalue  is  real 
and  distinct  and  it  can  easily  be  seen  that  all  the 
eigenvalues  are  less  than  1 and  Amax  approaches  1 
for  large  n . Thus  we  set  Amax  = 1 - e and  expand  Eq. 
(8)  to  lowest  order  in  e and  /',  which  gives  us 

+ (9) 

Then,  we  obtain 

Ca/l)  = A^ax  = (l-p/'"0",  (10) 

which  confirms  the  result  Eq.  (2)  by  the  sub- 
stochastic transition-matix  method. 

In  order  to  find  the  failure  probability  as  a func- 
tion of  applied  stress,  we  use  the  load-sharing  rule 
Eq.  (1),  coupled  with  the  fact  that  the  failure  of  the 
bond  carrying  the  largest  local  stress  nucleates 
catastrophic  failure,  and  thereby  use  the  relation 

oy,/(T=l-b|,  (11) 

where  at  is  the  breaking  strength  of  a single  fiber. 
Note  that  we  could  have  used  a variety  of  other 
load-sharing  rules  here,  and  for  example  the  same 
expression  with  n raised  to  an  arbitrary  power  is 
also  of  physical  interest.  This  result  combined  with 
Eq.  (10)  yields  the  probability  5(cr)  that  a fiber 
bundle  will  survive  at  stress  a 

Sia)  = (\-pf^-'y.  (12) 

For  large  n and  L,  this  becomes  the  modified 
Gumbel  form,  introduced  previously  [10,11]  in  the 
analysis  of  the  random  fuse  network.  Although 
CL{n)  in  Eq.  (10)  becomes  a Gumbel  distribution 
in  n,  the  substitution  of  n{a)  from  Eq.  (11)  pro- 
duces a modified  Gumbel  distribution  that  is  signif- 
icantly different  from  a Gumbel  form  in  cr  in  the 
high-reliability  tail  of  the  distribution.  This  modifi- 
cation is  discussed  further  in  Sec.  4. 

3.  Double-Cluster  Calculation 

Several  authors  [12,16,19,20]  have  suggested  that 
the  most  critical  defect  is  not  a single  cluster  of  n 
vacancies,  but  rather  a double  cluster  (double  co- 
linear  crack)  of  n vacancies  separated  by  a single 
occupied  site  located  anywhere  within  the  n + \ 
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adjacent  sites.  Such  a double  crack  is  shown  in  Fig. 
3a.  This  candidate  for  the  most  critical  crack  is  ap- 
pealing because  the  stress  enhancement  at  the  in- 
terior occupied  site  grows  as  n in  network  models 
rather  than  as  (as  for  the  edges  of  a single 
crack  in  a two-dimensional  network)  and  because 
the  increased  entropy  of  the  (n  + 1)  locations  of  the 
occupied  site  makes  it  more  probable.  Thus,  fol- 
lowing Harlow  and  Phoenix  [19],  we  consider  the 
probability  of  bundles  of  length  L with  repeated 
double  cracks  (and  single  cracks  when  the  occupied 
site  is  at  either  end)  not  exceeding  n vacant  sites  in 
any  two  adjacent  cracks  separated  by  a single  site. 
These  repeated  double  cracks  are  shown  in  Fig.  3b. 
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survival  configurations  (we  don’t  have  to  consider 
configurations  which  start  with  O’s). 

A great  simplification  in  the  characteristic  equa- 
tion 


<3) 
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X 


b) 


V. 


V 


j 


Fig.  3.  a)  A double  cluster  of  total  vacancies  plus  one  iso- 
lated occupied  bond,  b)  Repeated,  overlapping  double  clusters; 
each  pair  of  clusters  as  indicated  by  the  brackets  contains  € n 
total  vacancies  plus  one  isolated  bond. 


Harlow  [18]  showed  that  this  problem  is 
amenable  to  analysis  by  the  sub-stochastic  transi- 
tion-matrix method.  The  load-sharing  rule  is  still 
given  by  Eq.  (1)  as  before  but  now  n is  the  sum  of 
the  number  of  vacant  sites  immediately  on  the  left 
and  right  of  any  isolated  intact  bond  or  fiber.  Thus 
it  is  necessary  to  keep  track  of  not  only  the  number 
of  vacancies  in  the  cluster  being  considered  but 
also  those  in  the  previous  vacant  cluster.  There  are 
now  {n  +l)(n  +2)/2  distinct  endings  that  must  be 
considered  at  a site  (or  bundle  ending);  these  are 
given  by  the  basis  vector  <^  = <^(i),(/)(io)...(/)(io...o); 
<^(101),(A(101())-..<A(1()10...0);<^(1()01)...</>(!()010...0);-..;<^(10...01)’  where 
there  are  no  more  than  n total  vacancies  in  any 
element.  With  this  ordering  of  states  the  n =4  sub- 
stochastic matrix  for  this  problem,  for  example,  is 
given  by  Eq.  (13). 

With  periodic  boundary  conditions  Eq.  (5)  still 
holds  and  we  again  analyse  the  largest  eigenvalue 
of  A/„.  As  a technical  point,  note  that  since  we  are 
using  periodic  boundary  conditions,  we  can  always 
start  the  matrix  process  at  a surviving  bond,  and  so 
the  endings  considered  above  include  all  possible 


M„<f>  = \(f>. 


(14) 


where  <f>  are  the  eigenvectors  of  M„,  is  possible 
since  most  of  the  rows  of  M„  contain  only  a single 
non-zero  element  /.  This  gives,  for  example 
/</)(!)  = A</)(  10)  or  b4>(i)  = {flk)4i(i)  = (j)(w).  By  such  rela- 
tions, we  can  eliminate  all  the  rows  of  except 
the  rows  with  /?’s  corresponding  to  the  reduced  ba- 
sis vector  </>  = <^(i),(A(ioi),<^(iooi),...,(^(i(x)...oi).  The  result- 
ing (n  -f- 1)  equations  give  an  (n  -I- 1)  x (/j  -I- 1) 
matrb(  A/',  which  satisfies  the  reduced  characteris- 
tics equation 
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where  a =plk,  and  b =//A.  This  M'  matrix  can  be 
considered  as  a new  transition-matrix  which  adds  a 
cluster  at  a time  rather  than  a bond  or  fiber  at  a 
time.  Thus  we  have  the  characteristic  determinant 
equation. 
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For  small  n these  determinants  can  be  evaluated 
directly.  For  example, 

Dq  = — \-¥a,Di  = \—a  —a^b  and 


D2= -\+a+ab+aV-a^b\  (17) 

But  by  expanding  the  determinants  Eq.  (16)  by 
rows  and  columns,  we  can  show  that  there  is  an 
exact  recursion  relation 


bD„  (a,b)=  sDn  -2(ab,b)-  Dn-^{ab  \b),  (18) 

where  Dm{ab^,b)  is  Dn,{a,b)  with  a replaced  by 
{ab^,  and  where 

s = \+b-a^b''^\  (19) 

Note  that  s =s{a,b,n)=s{ab,b,n -2)=s{ab^,b,n 
— 4)  which  is  key  in  the  solvability  of  the  recursion 
relation  Eq.  (18).  After  some  detailed  analysis,  we 
have  found  (see  [21]  for  details),  that  this  recur- 
rence equation  may  be  solved.  The  resulting  char- 
acteristic equation  is  given  by 

7 (n-m 

, (20a) 


for  n ^3  and  odd.  While  for  n even  ^4,  we  find, 

nl2-l_^  na-l  2 
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In  Eqs.  (20a)  and  (20b), 
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The  key  quantity  s is  given  by  Eq.  (19)  above. 
Equations  (20a,b)  are  the  exact  expressions  for  the 
characteristic  equation  of  the  original  M in  Eqs. 
(13)  and  (14). 

Again  we  find  that  the  largest  eigenvalue  of  M is 
near  1 for  large  n , so  we  set  A = 1 — e and  expand 
Eqs.  (20a,b)  and  find  that  in  both  cases,  to  lowest 
order  in  e and  /”, 


from  the  possible  locations  of  the  single  bond  in  a 
double  cluster  of  size  {n  + 1).  The  ( —p ) is  a correc- 
tion to  properly  handle  the  single-cluster  cases  as 
well  as  the  double-cluster  case,  since  these  are  in- 
cluded whenever  the  isolated  bond  is  located  at 
one  end  of  the  double  cluster.  It  is  only  important 
for  small  n. 

In  order  to  check  and  better  understand  the 
asymptotic  form  Eq.  (21)  of  An,™  and  the  impor- 
tance of  the  other  eigenvalues  A,  we  have  made  sev- 
eral numerical  evaluations  of  the  various 
equations.  First,  we  have  numerically  found  the 
largest  eigenvalue  Amax  of  the  original,  full,  sub- 
stochastic transition-matrix  M as  given  by  Eq.  (13). 
Using  the  unit  vector  as  a starting  vector  we  re- 
peatedly apply  the  matrix  M to  it.  Since  the  largest 
eigenvalue  is  unique,  this  process  converges  expo- 
nentially to  the  largest  eigenvalue.  We  found  in 
general  that  convergence  occurred  to  six  significant 
figures  with  at  most  50  matrix  products  (even  for 
matrices  M of  dimension  {n  +l)(n  -I- 2)/2  = 10,000. 
The  sparsity  of  M with  this  iterative  procedure 
eliminated  matrix-storage  problems.  The  results  of 
this  iterative  procedure  for  Cz.(n)  versus  n are 
shown  in  Fig.  4 as  solid  lines,  and  the  dots  give  Amax 
with  Amax  as  given  by  the  asymptotic  form  Eq.  (21). 
Good  agreement  is  seen  for  all  p,  with  a small  devi- 
ation in  the  p-0.2  data  However,  a more  strin- 
gent test  is  needed  in  the  high-reliability  (large  n) 
tail.  Thus,  in  Fig.  5 we  plot  the  quantity 


1-A„ 


+ 1 yVl 


. + !==[(«  +2)p^-pl 


(22) 


n+ 1 


€=[in+2)p^-pf^^  + 0(f^’‘%  (21) 

Comparing  this  double-cluster  result  to  the  single- 
cluster result  Eq.  (9)  we  find  the  expected  (n  -t-  2) 


Fig.  4.  A plot  of  the  double-cluster  probability  Ctin)  vs  {n  + 1), 
for  p =0.2,  0.5,  and  0.8.  The  dotted  line  is  the  asymptotic  form 
given  by  Eq.  (24).  The  solid  lines  are  found  from  evaluating  the 
largest  eigenvalue  numerically,  and  by  using  C/,(n)  = Amax. 


341 


Volume  99,  Number  4,  July-August  1994 

Journal  of  Research  of  the  National  Institute  of  Standards  and  Technology 


Fig.  5.  Plots  of  elf"*^  vs  « + 1,  for  a) p =0.5  and  b) p =0.2.  The 
solid  lines  are  the  asymptotic  form  as  given  by  Eq.  (22);  the 
circles  are  the  exact  values  of  (1  - as  obtained  by  iter- 

ating M numerically. 

The  solid  straight  line  versus  (n  +1)  is  the  asymp- 
totic result,  which  is  linear  in  n and  this  is  com- 
pared with  the  iterated  numerical  values  (circular 
symbols),  for  /7=0.5  and  0.2.  For  large  n,  in  all 
cases  the  two  calculations  agree.  But  for  small  val- 
ues of  /j<0.5  there  appears  a minimum  in 
(1  - Amax)/T^‘  versus  n which  corresponds  to  higher 
order  terms  in/".  In  particular  the  next  order  term 
in  Amax  is  0(f^"'^)  which  would  appear  as  a 
correction  in  Fig.  5. 

Finally,  we  test  for  the  accuracy  of  neglecting  all 
but  the  maximum  eigenvalue  Amax.  It  is  possible  to 
directly,  numerically  evaluate  the  trace  of  by 
iteration  since  it  only  requires  storage  of  the  matrix 
and  one  vector  at  any  time.  Using  this  method,  we 
have  evaluated  the  quantity 

(23) 

which  should  converge  to  (1  - A^ax)//"^'  when  the 


largest  eigenvalue  is  dominant.  A numerical  test  of 
this  convergence  is  shown  in  Fig.  6,  for  L = 1000, 
and  shows  that  for  large  lattice  sizes,  the  most  im- 
portant corrections  to  Ctin)  are  the  higher  order 
contributions  to  Amax,  which  are  of  rather 

than  the  neglected  smaller  eigenvalues  of  M,  which 
are  relatively  unimportant  here. 

Finally,  we  obtain  the  asymptotic  form 

Cl  (n )-{!-[(«  +2)p'-/?r^'  + C(f'"'')}C  (24) 

Thus,  following  the  same  arguments  as  at  the  end 
of  Sec.  2,  we  find  that  to  leading  order  the  survival 
probability  of  the  entire  network  or  chain-of-bun- 
dles  is 

= (25) 


which  is  again  of  the  form  of  a modified  Gumbel 
distribution  with  slightly  different  coefficients  from 
Eq.  (12). 


Fig.  6.  Plots  of  the  exact  value  of  vs 

{n  + 1)  obtained  numerically  from  M (circles);  for  L = 1000,  for 
a)  p =0.5  and  b)  p =0.2.  The  solid  line  is  the  asymptotic  form 
Eq.  (22). 
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4.  Extreme-Statistical  Form 

For  large  L,  we  can  easily  find  the  limiting  form 
of  Clin)  as  given  by  Eqs.  (10)  and  (24)  respec- 
tively, 

Cr.(n)^~^exp[-Lp  exp(-n  log(l//))],  (26) 

and 

Clin)  ~ exp[-L((n  +2)p^ -p)exp{-n  log(l//))], 

(27) 

as  the  upper  limit  behavior  for  large  n which  is  im- 
portant in  the  high-surviveability  tail  of  the  distri- 
bution. This  is  a Gumbel  distribution,  as  is 
expected  from  the  exponential  [pf  or  npf]  be- 
havior for  the  probability  of  large  clusters  [Castillo, 
1988]. 

On  the  other  hand,  the  surviveabilities  S{a),  and 
5 (o-),  in  the  limit  L^oo,  are  found  to  be  modified 
Gumbel  forms  [from  Eqs.  (12)  and  (25)], 

SLia)~jxp[-Lp  exp(-(^-l)log(l//))]  (28a) 

and. 


ification  of  the  Gumbel  form  is  essential  in  order  to 
correctly  represent  the  important  high-reliability 
tail.  Note  that  this  is  not  true  of  the  Weibull  distri- 
bution, which  as  well  as  being  a stable  limiting  ex- 
treme-value distribution,  does  have  the  physically 
correct  behavior  as  stress  approaches  zero.  This  is 
one  good  reason  why  the  Weibull  distribution  is  a 
very  robust  form  in  the  analysis  of  failure  problems. 
We  suggest  that  the  family  of  modified  Gumbel 
distributions  of  the  sort  Eq.  (28),  should  be  simi- 
larly robust,  in  contrast  to  the  conventional  Gum- 
bel distribution  which  is  of  limited  use  in  the 
analysis  of  the  statistics  of  material  failure. 

5.  General  Scaling  Behavior 

The  size  dependence  and  general  form  of  the 
limiting  distribution  can  usually  be  found  from  a 
back-of-the-envelope  scaling  calculation  which  we 
introduced  previously  [9,11].  First,  for  the  single- 
cluster calculation,  consider  the  probability  Pein) 
of  finding  a cluster  of  size  n in  a sample  of  size  L. 
The  order  of  magnitude  of  this  probability  is 

PUn)  « Lp^-r/p=Lpr,  (32) 


SL(cr)^~jxp[-Li^p^-p ) exp(  - (^  - 1) 

log(l//))].  (28b) 

That  the  dominant  behavior  of  5z.(cr),  and  Sl{o-)  is 
exp(—LA  exp(— f))  as  cr  tends  to  zero  is  essential 
to  ensuring  that  the  survival  distributions  have  the 
proper  limiting  approach  to  one  when  the  applied 
stress  approaches  zero.  Flarlow  and  Phoenix  [19] 
have  numerically  shown  that  this  high-reliability 
tail  can  not  be  well  described  by  a Gumbel  form  for 
a (such  as  exp( exp(5((T  - cro)))).  But  this 
failure  is  obvious  since  the  Gumbel  form  doesn’t 
approach  one  until  o-^—  oo,  so  at  sufficiently  small 
stresses  it  must  be  inaccurate.  Nevertheless  the 
standard  texts  on  extreme  distributions  (see,  for  ex- 
ample [8]  or  [22])  seem  to  suggest  that  the  Gumbel 
distribution  is  the  appropriate  one  in  such  cases. 
The  difference  clearly  is  in  the  form  of  the  nor- 
mally assumed  scaling  limAf-,*[5 (o-)T  = S {aNO-  + bN), 
which  fits  the  shift  and  slope  of  the  limiting  func- 
tion S{a)  at  its  median,  but  fails  near  the  high-reli- 
ability limit  cr  = 0.  It  would  seem  that  instead  a 
more  general  scaling  form  limAf-^=c[S(n(cr))T  = 
5(fl/vn(cr)-t-6/v)  must  be  allowed  to  also  include  the 
proper  high-surviveability  limit  near  zero  stress.  In 
many  practical  material-failure  problems,  this  mod- 


since, for  normalization,  Xp^f'—p,  and  since 
there  are  L different  locations  in  the  sample  where 
a cluster  could  be  located.  For  the  maximum  clus- 
ter size  to  be  expected  in  a sample  of  size  L we  set 

P,(n)~Lpf--^l,  (33) 

and  obtain  the  size  dependence 


''"^*-log(V/)’ 


(34) 


or,  from  the  load-sharing  rule  Eq.  (1),  the  break- 
down stress  Ob 


’ (35) 

1 "1"  ^max/2  ^ . In  L 

21n(l//) 

scales  to  zero  logarithmically  in  the  thermodynamic 
limit.  A similar  argument  for  double  clusters  gives 

PL{n)<^L{n+\)p^f\  (36) 

since  there  are  (n+1)  places  to  put  an  isolated 
bond  in  the  /i -double-clusters.  We  then  obtain  the 
same  limiting  form  Eq.  (35)  for  a-c(L),  although 
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there  are  additive  (log(logL))  corrections  in  the 
double-cluster  case.  The  logarithmic  scaling  law  in 
turn  implies  that  the  failure  statistics  is  of  the  dou- 
ble-exponential form  given  in  Eq.  (28).  The 
Weibull  and  Frechet  distributions  always  give 
power-law  size  scaling.  These  qualitative  argu- 
ments are  very  powerful  and  are  confirmed  by  the 
rather  elaborate,  exact  calculation  described  here. 
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